ALGEBRO-GEOMETRIC SOLUTIONS OF A DISCRETE SYSTEM 
RELATED TO THE TRIGONOMETRIC MOMENT PROBLEM 



JEFFREY S. GERONIMO, FRITZ GESZTESY, AND HELGE HOLDEN 



Abstract. We derive theta function representations of algebro-geometric so- 
lutions of a discrete system governed by a transfer matrix associated with (an 
extension of) the trigonometric moment problem studied by Szego and Bax- 
ter. We also derive a new hierarchy of coupled nonlinear difference equations 
satisfied by these algebro-geometric solutions. 



1. Introduction 

Let {a(n)}„gN C C be a sequence of complex numbers subject to the condition 

|Q!(n)| < 1 for all neN, (1.1) 
and define the transfer matrix 

r(.).(4 f), ..T, (1.2) 

with spectral parameter z on the unit circle T = {zgC||z| = 1}. Consider the 
system of difference equations 

$(z,n) ^ T(z,n)$(z,n- 1), (z,n)eTxN (1.3) 

with initial condition $(z, 0) = ( } ), z e T, where 

H^, n) = > (z, n) e T X Nq. (1.4) 

(Here Nq = N U {0}.) Then <y9( • , ti) are monic polynomials of degree n and 

(p*(z,?i) = z"^(l/z,n), (z,n)eTxNo, (1.5) 

the reversed *-polynomial of f{z, n), is of degree at most n. These polynomials were 
first introduced by Szego in the 1920's in his work on the asymptotic distribution 
of eigenvalues of sections of Toeplitz forms [40], [41] (see also [32, Chs. 1-4], [42, 
Ch. XI]). Szego's point of departure was the trigonometric moment problem and 
hence the theory of orthogonal polynomials on the unit circle: Given a probability 
measure da supported on an infinite set on the unit circle, find monic polynomials 
of degree n in z = e*^, 6 E [0, 27r], such that 

r2iT 

j{n fda{e'^) <^{e^\ m)^{e'^,n) = <5™,„, m, n € No, (1.6) 





Date: July 23, 2004. 

1991 Mathematics Subject Classification. Primary 35Q53, 58F07; Secondary 35Q51. 
The research of the second and third author was supported in part by the Research Council of 
Norway. 

Supported in part by the US National Science Foundation under Grants No. DMS-0200219 
and DMS-0405526. 



2 



J. GERONIMO, F. GESZTESY, AND H. HOLDEN 



where 

/ x2 1 1 for n = 0, ,^ . 

[rij-i (1- |a(j)r) forneN. 
Here we chose to emphasize monic polynomials ip{- ,n) in order to keep the factor 
7 out of the transfer matrix T. Szego showed that the polynomials (1.4) satisfy 
the recurrence formula (1.3). Early work in this area includes contributions by 
Akhiezer [9, Ch. 5], Geronimus [24], [25], [26, Ch. I], Krcin [33], and Tomcuk 
[43]. For a modern treatment of the theory of orthogonal polynomials on the unit 
circle and an exhaustive bibliography on the subject we refer to the forthcoming 
monumental two- volume treatise by Simon [38] (sec also [39]). For fascinating 
connections between orthogonal polynomials and random matrix theory we refer, 
for instance, to Dcift [18]. 

An extension of (1.3) was developed by Baxter in a series of papers on Toeplitz 
forms [10]-[13] in 1960-63. In these papers the transfer matrix T in (1.2) is replaced 
by the more general transfer matrix 

^(^)-(;, tj (1.8) 

with a = a(n), (3 = P{n), subject to the condition 

a{n)(3{n) ^ 1 for all n G N. (1.9) 

Studying the following extension of (1.3), 

^{z, n) = U{z, n)^{z, n - 1), {z, n) G T x N, (1.10) 

Baxter was led to biorthogonal polynomials on the unit circle with respect to a 
complex-valued measure. In this paper we will primarily be concerned with Baxter's 
extension (1.10) of (1.3). 

To simplify our notation in the following, shifts on the lattice are denoted using 
superscripts, that is, we write for complex-valued sequences /, 

(5±/)(n) = /±(n) = f{n ±1), n e Z, {/(n)}„ez C C (1.11) 

and apply the analogous convention to 2 x 2 matrices and their entries. 

In the mid seventies, Ablowitz and Ladik, in a series of papers [3], [4], [5], [6] 
(see also [1], [2, Sect. 3.2.2], [7, Ch. 3]), used inverse scattering methods to analyze 
certain integrable differential-difference systems. One of their intcgrable variants 
of such systems, a discretization of the AKNS-ZS system, is of the type 

-iat - (a+ -2a + a") + aP{a+ + a") = 0, (1.12) 

-iPt + (/?+ - 2/3 + /?-) - aP{(3+ +/?-)= (1.13) 

with a = a{n, t), f3 = (3{n, t). In particular, Ablowitz and Ladik [4] (see also [7, Ch. 
3]) showed that in the focusing case, where (3 = —a, and in the defocusing case, 
where (3 = a (cf. (1.2)), (1.12) and (1.13) yield the discrete analog of the nonlinear 
Schrodinger equation 

-iat + 2q! - (1 ± |a|^)(Q;+ + a") = 0. (1.14) 

Algebro-geometric solutions of the AL system (1.12), (1.13) have been studied 

by Ahmad and Chowdury [8], Bogolyubov, Prikarpatskii, and Samoilenko [15], 
Bogolyubov and Prikarpatskii [16], Ceng, Dai, and Cewen [20], Vekslerchik [44], 
and especially, by Miller, Ercolani, Krichever, and Levermore [34] in an effort to 
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analyze models describing oscillations in non-linear dispersive wave systems. In [34] 
the authors use the fact that the AL system (1.12), (1.13) arises as the compatibility 
requirements of the equations 

$ = [/$~, $t- = VK$-. (1.15) 

Here U is precisely Baxter's matrix in (1.8) and W is defined as follows, 

= (;, °) , ^(.) = . {' , I - °t:_.) . (1.16) 

Thus, the AL system (1.12), (1.13) is equivalent to the zero-curvature equations 

Ut + UW -W+U = Q. (1.17) 

Miller, Ercolani, Krichever, and Lcvermorc [34] then performed a thorough analysis 
of the solutions $ = $(z,n, t) associated with the pair {U,W) and derived the 
theta function representations of a,(i satisfying the AL system (1.12), (1.13). In 
the particular focusing and defocusing cases they also discussed periodic and quasi- 
periodic solutions a with respect to n and t. 

Unaware of the paper [34], Geronimo and Johnson [22] studied the defocusing 
case (1.3) in the case where the coefficients a are random variables. They provide a 
detailed study of the corresponding Weyl-Titchmarsh functions, m± , which satisfy 
the Riccati-type equation (for z G C\T, n e Z), 

a{n)m±{z, n)m±{z, n — 1) — m±{z,n— 1) + zm±{z, n) = za{n) (1-18) 

(which should be compared to the identical equation (3.20) for the fundamental 
function </> in the defocusing case (i = a). These functions take on the values 
\m+{z)\ < 1 and |m_(^)| > 1 for [^l < 1 (cf. [21], [22]). Utilizing the fact that m+ 
is a Schur fimction (i.e., analytic in the open imit disc with modulus less than one) 
and the close relation between such functions and the orthogonality measure da+, 
they perform the transformation 

$ = A$, A=^r. ~.M . (1.19) 



^/2 \i i 

With this change of variables ra± transform into 

m±(^,n) =i^i^^^4^. 0eC\T, neZ. (1.20) 
1 — m±{z,n) 

The Schur property of m+ (equivalently, the relation between Schur functions, 
Caratheodory functions, and positive measures on the unit circle [9], [37]-[39]) 
implies the standard representation, 

fh+{z, n) = i da+(e'^ n) , z € C\T, n € Z. (1.21) 

Jo — z 

Under appropriate crgodicity assumptions on a and the hypothesis of a vanishing 
Lyapunov exponent on the prescribed spectral arcs on the unit circle T, Geronimo 
and Johnson [22] showed that the m-functions associated with (1.3) are reflection- 
less, that is, ?7i_|- is the analytic continuation of through the spectral arcs and 
vice versa, or equivalently, rh± are the two branches of an analytic function m on 
the hyperelliptic Riemann surface with branch points given by the end points of the 
spectral arcs on T. They developed the corresponding spectral theory associated 
with (1.3) and the unitary operator it generates in £^(Z) (cf. [23]). This can be 
viewed as analogous to the case of real-valued finite-gap potentials for Schrodinger 
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operators on R (cf., e.g., [14], [27]) and self-adjoint Jacobi operators on Z (of., e.g., 
[17]). In particular, Geronimo and Johnson [22] prove the quasi-periodicity of the 
coefficients a in the defocusing case /3 = S. Connections with aspects of integrabil- 
ity, a zero-curvature or Lax formalism, and the theta function representation of a, 
are not discussed in [22] . The whole topic has been reconsidered in great detail and 
partially simplified in the upcoming two- volume monograph by Simon [38, Ch. 11] 
and aspects of integrability (Lax pairs, etc.) in the periodic defocusing case will 
further be explored by Nenciu and Simon [36]. 

The principal contribution of this paper to this circle of ideas is a short derivation 
of theta function formulas for algcbro-geometric coefficients a, /3 associated with 
Baxter's finite difference system (1.10). Rather than considering solutions of a 
particular AL flow such as(1.12), (1.13), we will focus on a derivation of the coupled 
system of nonlinear difference equations satisfied by algebro-geometric solutions 
Q!,/3 of (1.10) (a new result) and its algebro-geometric solutions. In this sense 
our contribution represents the analog of determining algebro-geometric coefficients 
(generally, complex-valued) in one-dimensional Schrodinger and Jacobi operators 
and deriving the corresponding Its-Matveev-type theta function formulas. As a 
by-product in the special defocusing case (3 = a with |a(n)| < 1, n G Z, we recover 
the original result of Geronimo and Johnson [22] that a is quasi-periodic without 
the use of Fay's generalized Jacobi variety, double covers, etc. 

In Section 2 we describe our zero-curvature formalism and the ensuing hierarchy 
of nonlinear difference equations for a, /3. Our principal Section 3 then is devoted 
to a detailed derivation of the theta function formulas of all algebro-geometric 
quantities involved. Appendix A collects relevant material on hyperelliptic curves 
and their theta functions and introduces the terminology freely used in Section 3. 

2. Zero- Curvature Equations and Hyperelliptic Curves 

In this section we introduce the basic zero-curvature setup for algebro-geometric 
solutions of (1.10). We follow the approach employed in [17], [27]-[30] in the anal- 
ogous cases of stationary KdV, AKNS, and Toda solutions. 

We start by introducing the complex-valued sequences 

{a(n)}„6z,{/3(n)}„ezCC, (2.1) 

and define the recursion relations 

/o = -2a+, 50 = 1, ho = 2(3, (2.2) 

gi+i - 57+1 = ahj + pfe, £ e No, (2.3) 

fi+i = ft- a{ge+i + 57+1), i e No, (2.4) 

he+i = hj + I3{ge+i + 57+1), £ e No, (2.5) 

Here shifts on the lattice arc denoted using superscripts as introduced in (1.11). 
In addition we get the relations 

- 9i+i = ahg+i + Pff+i, £ G No, (2.6) 

which are derived as follows, 

ahe+i + I3fl^i = ahj + aP{ge+i + gj^-^) + (3fe - af3{g£+i + J 

= ahj +pfe = ge+i-gj^+i, ^ e Nq, (2.7) 

using relations (2.4), (2.5), and (2.3). 
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Remark 2.1. One can compute the sequences {fg}, {gi}, and {he} recursively 
as follows. Assume that fe, gi, and are known. Equation (2.3) is a first order 
difference equation in gc+i that can be solved directly and yields a local lattice 
function. The coefficient ge+i is determined up to a new constant denoted by 
Q+i e C. Relations (2.4) and (2.5) then determine fe+i and he+i, etc. 

Explicitly, one obtains 

/o = -2a+, /i = 2{{a+fp + a+a++(3+ - a++) + ci(-2a+), 

go = l, gi = -2a+l3 + ci, (2.8) 

ho = 2(3, hi = 2{- - aP'P + /?") + ci2/3, etc., 

where {ce}e^jq C C denote certain summation constants. 

Next, assuming ^ e C, we introduce the 2x2 matrix U (z) by 

U{z,n)=( J", , a{n)\ ^ ^ ^_ ^^.g) 



^zP{n) 1 

In addition, we introduce for each fixed p G N the following 2x2 matrix V^_|_i(z), 

supposing Fp{- ,n) and Gp+i{- ,n), i?p+i(- ,n) to be polynomials of degree p and 
p+l, respectively (cf., however. Remark 3.1), with respect to the spectral parameter 

Postulating the stationary zero-curvature condition 

U{z,n)Vp+i{z,n)-Vp\,{z,n)U{z,n) = 0, p e No, (2.11) 

then yields the following fundamental relationships between the polynomials Fp, 
Gp+i, and Hp+i, 

Fp-zFp -a{Gp+i + G;+,)=0, (2.12) 
zl3{Gp+i + G;+i) + H-^, - zHp+i = 0, (2.13) 
z{G-^^ - Gp+i) + aH-^^ + zf3Fp = 0, (2.14) 
Gp+i - - aHp+i - zpFp = 0. (2.15) 

Moreover, using relations (2.12)-(2.15) one shows that the quantity Gp+i —FpHp^i 
is a lattice constant and hence the expression 

G'p+i(z,n)^ - Fp{z,n)Hp+i{z,n) = i?2p+2(z) (2.16) 

is an n-independent polynomial of degree 2p-\-2 with respect to z. (That Gp+i — 
FpHp+i, 7^ 0, is a lattice constant also immediately follows from (2.11) taking 
determinants.) 

In order to make the connection between the zero-curvature formalism and the 
recursion relation (2.2)-(2.5), we now introduce the polynomial ansatz with respect 
to the spectral parameter z, 

p p+l p+l 

Pp{^) = ^fp-iz\ Gp+i{z) = ^gp+i-iz\ Hp+i{z) = ^hp+i-iz^. (2.17) 
t=0 i=0 c=o 
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The stationary zero-curvature condition (2.11) imposes further restrictions on the 
coefficients of Vp+i that we will now explore. Since 170 = 1, the quantity i?2p+2 in 
(2.16) is a monic polynomial of degree 2p + 2, that is, 

2p+l 

R2p+2{z) = Y\_i^~-^m), {Em}m=0,...,2p+1 C C (2-18) 
m=0 

Next we assume p € N to avoid cumbersome case distinctions concerning the trivial 
case p = 0. Insertion of (2.17) into (2.12)-(2.15) then yields the relations (2.2) 
(normalizing go = 1) and the recursion relations (2.3), (2.4), and (2.5) for ^ = 
0, . . . ,p — 1. In addition, one obtains the equations 



fp - a(.9p+i 4 


-9p+i) 


= 0, 


(2.19) 


0i9p+i + 9p+i) + hp - 


- hp+i 


= 0, 


(2.20) 




hp+i 


= 0, 


(2.21) 


Qp+i - gp+i + ahp 


+ Pfp 


= 0, 


(2.22) 


ahp+i +9p+i- 


- 9p+i 


= 0. 


(2.23) 



Moreover, one infers the relations (cf. (2.6)) 

gi-gi =ahe + pfi, £ = 0,...,p. (2.24) 

Combining (2.21) and (2.23), we first conclude that gp+i is a lattice constant, 
that is, 

gp+i = gp+i e c. (2.25) 

In addition, using (2.20), (2.21), and (2.25) one obtains 

= hp+i =h-+ (5{gp+i + 5-+1) = h-+ 2gp+iP, (2.26) 

and hence, hp = — 2,gp4.i/3+. Equations (2.19) and (2.25) also yield fp = 2gp+ia in 
agreement with (2.22). Moreover, (2.25) is consistent with taking 2; = in (2.16) 
which yields 

2p+l 

9I+1 = n (2.27) 

Thus, the stationary zero-curvature condition (2.11) is equivalent to a coupled 
system of nonlinear difi'erence equations for a and /3 which we write as 

s-SBpia,P) = = 0, gp,. = g^,,, (2.28) 

in honor of the pioneering work by Szcgo and Baxter in connection with the transfer 
matrices (1.2) and (1.8). Varying p G No in (2.28) then defines the corresponding 
stationary SB hierarchy of nonlinear difference equations. The first few equations 
explicitly read 

s-SBo (a, (3) = 2 (^"".^^"^^Jj") =0, 51 = 

s-SBi(a,/3) - 2 [_^-f3-f3- _ ^(^-)2 + - 0, (2.29) 

g2 = gt, etc. 
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By definition, the set of solutions of (2.28), with p ranging in No, represents the 
class of algcbro-gcometric solutions associated with Baxter's finite difference system 
(1.10). The hierarchy of coupled nonlinear difference equations (2.28) is new. 

Remark 2.2. [i) The scaling behavior fe.{Aa, A"\/3) = Af/Xa, /?), g/XAa, A-'^(3) = 
ge{a,P), he{Aa,A-^(3) = A-^he{a,l3), leNo, Ae C\{0}, shows that the station- 
ary SB hierarchy (2.28) has the scaling invariance, 

{a, p) ^ {Aa, A-^13), A e C\{0}. (2.30) 

In the special focusing and dcfocusing cases, where (3 = —a and P = a, respectively, 
the scaling constant A in (2.30) is further restricted to 

1^1 = 1. (2.31) 

(m) In the defocusing case (3 — a, the compatibility requirement of the two equa- 
tions in (2.28) requires the constraint Iffp+iP = 1 and additional spectral theoretic 
considerations in connection with the trigonometric moment problem, assuming 
\a{n)\ < 1, n G Z, enforce {E„i}m=o,...,2p+i C T. The additional condition of 
periodicity of a then implies further constraints on {Ejn}m=o,...,2p+i (cf. [38, Ch. 
11]). The special case of real-valuedness of a also enforces additional constraints 

on {Em}m=0,...,2p+1- 

3. Theta Function Representations 

In this our principal section, we present a detailed study of algebro-geometric 
solutions associated with (1.10) with special emphasis on theta function represen- 
tations of a, (3 and related quantities. Wc employ the techniques discussed in [17] 
and [27] in connection with other integrable systems such as the KdV, AKNS, and 
Toda hierarchies. 

Throughout this section we suppose 

{a(n)}„ez, {/3(n)}„ez C C, a{n)P{n) ^ 0, 1, n G Z, (3.1) 

and assume (2.2)-(2.5), (2.11), (2.17). Moreover, we freely employ the formalism 
developed in Section 2, keeping p e No fixed. 

Returning to (2.18) we now introduce the hyperelliptic curve Kp with nonsingular 
affine part defined by 

tZp : Tp{z, y) = y^ - R^p+^iz) = 0, (3.2) 

2p+l 

i?2p+2(-2) = n - ^™)' {'^m}m=0....,2p+l C C\{0}, (3.3) 

m— 

Em 7^ Era' for m ^ m', m, m' = 0, . . . , 2p + 1. (3-4) 

Equations (3.1)-(3.4) are assumed for the remainder of this section. We compactify 
/Cp by adding two points Poo+ and Poo_ , -Poo+ 7^ -Foo_ , at infinity, still denoting its 
projective closure by Kp. Finite points P on Kp are denoted by P = [z^y) where 
y{P) denotes the meromorphic function on Kp satisfying Tp{z, y) = 0. The complex 
structure on ICp is then defined in a standard manner and )Cp has topological genus 
p. Moreover, we use the involution 

*:/Cp^/Cp, P= (z,y)^P* = (z,-?y), Poc± ^Poc±* =Poo^. (3.5) 

For further properties and notation concerning hyperelliptic curves we refer to Ap- 
pendix A. 
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Remark 3.1. The assumption a{n) ^ 0, /3(n) 7^ 0, n e Z, in (3.1) is not an 
essential one. It has the advantage of guaranteeing that for all n S Z, Fp(- ,n) 
and iJp+i( • , n) are polynomials of degree p and p + 1, respectively. If a"'"(no) = 
(resp., /3(no) = 0) for some no € Z, then Fp( • , no) has at most degree p — 1 (resp., 
i?p+i(- ,no) has at most degree p). The latter n-dependence of the degree of the 
polynomials and iJp+i enforces numerous case distinctions in connection with 
our fundamental function in (3.14) below. For simplicity we will in almost all 
sitiiations avoid these cumbersome case distinctions and hence assume a(n) 7^ 0, 
/3(n) 7^ 0, n e Z throughout this section. (For an exception see Remark 3.5.) In the 
extreme case that a = (i.e., a(n) = for all n G Z), then Fp = and hence the 
curve /Cp becomes singular as ii2p+2(-2) = Gp+i(z,n)^, z G C, by (2.16), and thus 
the branch points of /Cp necessarily occur in pairs. (In addition, Gp+i(2;,n) then 
becomes independent of n G Z.) The same argument applies to /3 = since then 
i/p_l_i = 0. For this reason the trivial cases a = and /3 = in (3.1) are excluded in 
the remainder of this paper. Finally, in order to avoid numerous case distinctions 
in connection with the trivial case p = 0, we shall assume p > 1 for the remainder 
of this section (with the exception of Example 3.13). 

In the following, the zeros of the polynomials Fp( • , n) and i?p+i( • , n) (cf. (2.17)) 
will play a special role. We denote them by {A*j('^)}j=i,...,p and {t'^(n)}^=o,...,p and 
hence write 

p p 
Fp(z) = -2Q+[](^-M,), i?p+i(z) = 2/3[](2-i.,). (3.6) 

In addition, we lift these zeros to /Cp by introducing 

iij{n) = {iJ,j{n},Gp+i{iJ,j{n),n)) e /Cp, j = 1, . . . ,p, n e Z, (3.7) 
ue{n) = {ue{n),-Gp+i{i^i{n),n)) e ICp, £ = 0, . . . ,p, n € Z. (3.8) 

We recall that hp+i = (cf. (2.21)). Hence we may choose 

z/o(n) = 0, n e Z. (3.9) 

Define 

Po,± = (0, ±Gp+i(0)) = (0, ±5p+i), (3.10) 

where 

2p+l 

2/(Po,±) = ±ffp+i, 5p+i = n (3.11) 

m=0 

Wc emphasize that Po,± and Poo± arc not necessarily on the same sheet of /Cp. The 
actual sheet on which Po,± lie depends on the sign of gp+i- Thus, one obtains 

i>o = Po,-. (3.12) 

The branch of j/( • ) near Poo± is fixed according to 

lim Tr^^T^TT^= lim ^'i^L = Tl- (3.13) 
|.(p)Hoo Gp+i(z(P)) |.(P)Hoo z(P)P+i ^ ^ ' 

Next, we introduce the fundamental mcromorphic function </> on /Cp by 

^(P.n) = ^ + ^f^^!'"^ = P=(.,,)GX:p,nGZ (3.14) 

Pp(z,n) 2/-Gp+i(2;,n) 
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with divisor ((/>( • , n)) (cf. the notation for divisors introduced in (A. 20) and (A. 21)) 
given by 

• ,n)) = X>pj, _^(„) - 2?p„_^(n). (3.15) 
Here we abbreviated (cf. (A.20), (A.21)) 

/£ = {/ti, . . . ,/ip}, V = {j>i, ...,i>p}e SymP ICp. (3.16) 
The stationary Baker- Akhiezer vector \1/(P, n, no) is defined on ICp by 

*^^'"'"°^=U(^,n,no)j' (^-1^) 

fni^=„o+i(^ + "(H'/'"('P>'Ti)). n>no + l, 
■^i(P,n,no) = < 1, n = nQ, (3.18) 

inl=„+i(^ + a(m)r (P,"^))-\ n < no - 1, 

[ ni=„o+i(^/3(m)r (i^,m)-i + 1), n > no + 1, 

V'2(-P,n,no) = (?i'(-P,no) < 1, n = no, 

in;=„+i(^/5Mr (J'-m)-' + l)-\ n < no - 1, 

P e /Cp, (n,no) e (3.19) 

Clearly \E'(-,n, no) is mcroinorphic on )Cp since 0(-,rn) is meromorphic on /Cp. 
Fundamental properties of (p and ^' arc summarized next. 

Lemma 3.2. Suppose a, f3 d <C satisfy (3.1) and the pth stationary SB system 
(2.28). Moreover, assume (3.2)-(3.4) and let P = {z,y) G /Cp\{-Pcx)+)-Poo_}) 
(n, no) e I?. Then cf) satisfies the Riccati-type equation 

a(j){P)(j)- (P) - (l>- (P) + z(l>{P) = zp (3.20) 

and ^ fulfills 

MP,n,no)/MP,n.ri,o) = (l){P,n), (3.21) 

*(P,n,no) = i7(z,n)*-(P,n,no), (3.22) 

-y*-(P,n,no) =yp+i(z,n)*-(P,n,no). (3.23) 

Proof Using y^ = Gl+^ - FpHp+i (cf. (2.16), (3.2)) and (3.14), the left-hand side 
of (3.20) can be rewritten as follows 

a#- -cj)- +z(P-z(3^ {FpFp)-^ ["(Gp+i - FpHp+i + y{Gp+i + G^+J 

+ Gp+iGp+i) - (2/ + G-^^)Fp + z{y + Gp+,)F; - zpFpP-] . (3.24) 

Insertion of (2.12) and (2.15) into (3.24) then shows that the right-hand side of 
(3.24) vanishes. This proves (3.20). Relation (3.21) is proven inductively as follows. 
Since it holds for n = no by (3.18), (3.19) wc assume that 

i}2iP,m,no)/Tpi{P,m,no) ^ (j){P,m), m = uq, . . . ,n - 1. (3.25) 

Then combining (3.18), (3.19), and (3.25), one obtains 

MP,n,no) ^ zp{n)cp-{P,n)-' + l 

MP,n,no)'^ ^ '""^ z + ain)cf>-{P,n) ' ^""-^^^ 

and hence 

■!/'2(P,n,no) il)2{P,n,no) 

'^WtT^ {P,n)-(p {P,n)+z—-^ ':--zP{n) = 0. (3.27) 
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Comparison with (3.20) (cf. also (3.24)) then proves (3.21) for all n > uq. The case 
n < no — 1 is proven analogously. By (3.18) and (3.19) one infers 

V'l {P, n, no) = [z + a{n)(j)~ {P, n)]V'r {P, n, no) 

= z'4)'^{P,n,nQ) + Q;(n)V'^(P,n,no), (3.28) 

V'2(P,n,no) = [zl3{n)(l>-{P,n)-^ + l]V2"(f, n, no) 

= z/3(n)V^f (P, n, no) + ^J^{P, n, no), (3.29) 
by (3.21). This proves (3.22). An apphcation of (3.14) implies 

I'p+l* — \ rr- J,- _ n- ./,- ~ WW- ('^-^-l/^- ^ 



(3.30) 

and hence (3.23). □ 

We note that the Riccati-typc equation (3.20) for <p coincides with that of m± 
in (1.18) in the defocusing case (3 = a. 

Next, we derive trace formulas for a and /3 in terms of the zeros and Vj of Fp 
and Hp+i, respectively. For simplicity we just record the simplest case below. 

Lemma 3.3. Suppose that a,(3 cC satisfy (3.1) and the pth stationary SB system 
(2.28). Then, 

JL^tirifl,, ^^tirifi.. (3.31) 

a+ 9p+i fJi f3 9p+i fj^ 

Proof. Combining (2.17), fp = 2gp+ia, and (3.6) yields 

p p 
2gp+ra = fp = fo{-ir n = -2a+(-lF H l^j' (3-^2) 

Using hp = — 2£(p+i/3+, the case P~^/l3 is analogous. □ 

The following result describes the asymptotic behavior of (/> as P ^ Poo± and 
P ^ Po,±. 

Lemma 3.4. Suppose that a, /3 C C satisfy (3.1) and the pth SB system (2.28). In 
addition, assume (3.2) -(3.4) and let P = {z,y) G )Cp\{Poo+,Poo-} o-nd (n, no) € 
I?. Then, 



/3+(l-a/3)/3-C + 0(C') asP^Poo^, 



C-o |_(a+)-iC-^ + (l-a+/3+)a++(a+)-2 + 0(C) as P ^ Poo_, 

C = 1/2, (3.33) 

MP^ - -(l-«/3)«-WC + 0(C2) asP^Po,+ , 

'^^^^f-o|-/3+C-(l-a+/3+)/3++C^ + 0(C3) asP^Po,-, ^~ ^ 

Proof. Inserting the ansatz 

[Mn) + Mn)z + 0{z'^) asz^O ^' ' 

into the Riccati-type equation (3.20) produces (3.33) and (3.34). □ 
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Using (3.14)-(3.19) one obtains for the divisor {ipj {-,71, no)) of the meromorphic 
functions Vj( " > f^t '^o), j = 1, 2, 

(V'i(- ,n,no)) = ©^(„) - I'^(„o) + in- no){Vp^ _ - Vp^^), (3.36) 

(V2(-,n,no)) =Pp„,_^(„) -Pp^_^(„„) + (n-no)(DPo,_ -Pp^^). (3.37) 

Next, we briefly consider the asymptotic behavior of (j) in the case where the 
conditions a{n)l3{n) ^ are violated for some n G Z. 

Remark 3.5. First we note that if a+ ^ 0, then by (3.33) no pole fij of (p hits the 
point Poo_- Similarly, by (3.34), Po,_ = i>o is a zero of (j). The case /3(n) = for 
some n G Z poses no difficulty and (3.33) and (3.34) extend continuously in this 
case. The case a(n) = for some n e Z is more involved and causes higher order 
poles in (3.33) and (3.34). An expHcit calculation yields — 1/z) 

0(P) = l^^^l^^ ^ , asP^Poo+ if = 0, a++ ^ 0. (3.38) 

|-(a++)-iC-' + 0(C"') asP^Poo_ ^ y J 

Thus, if Q!+ = 0, 0"*"+ ^ 0, one of the poles fij of (j) hits the point . However, 
still no pole of (p hits Poo+ . Similarly, using 

fp = 25p+ia, fp-i = 2gp+i{a~ - o?l3+ - a~ap) + 2Ca, 

9p = -2gp+iaf3+ + g;l^2{2p + l)ci, (3.39) 

hp = -2gp+ip+, hp-i = 2gp+i{-l3++ + app+ + a- 0") - 2C/3, 

one derives {C, = z) 

HP) = + asP^Po,. if,^o,.-^0. (3.40) 

\-o|0(C) asP^Po,- ' ^ ^ ' 

Thus, if Q = 0, 7^ 0, one of the poles fij of (p hits the point Po,+. In addition, 
Pq,- remains a zero of (p. 

Since nonspecial divisors will play a fundamental role in this section, we now 
take a closer look at them. 

Lemma 3.6. Suppose that a, f3 cC satisfy (3.1) and the pth stationary SB system 
(2.28). Moreover, assume (3.2)-(3.4) and let n G 7,. Let Vp,, jj, = {/ti, . . . , /ip} 

and Dp, v_ = {i^i, . . . ,i>_p}, be the pole and zero divisors of degree p, respectively, 
associated with a, (3 and (p defined according to (3.7), (3.8), that is, 

fij{n) ^ {fij{n),Gp+i{fij{n),n)), j = 1, . . . ,p, n e Z, (3.41) 

Pji^n) = {iyj{n),-Gp+i{iyj{n),n)), j = 1, . . . ,p, n e Z. (3.42) 

Then 2?/i(„) and I?j>(„) are nonspecial for all n gZ. 

Proof. We provide a detailed proof in the case of By Theorem A. 2, 

is special if and only if {fii{n), . . . ,fip{n)} contains at least one pair of the type 
{/t(n), /i*(n)}. Hence f is certainly nonspecial as long as the projections /ij(n) 
of fij{n) are mutually distinct, /Uj(n) ^ /Ufe(n) for j ^ k. On the other hand, if two 
or more projections coincide for some no € Z, for instance, 

Wi (no) = ---=IJ'jN (no) = fJ-o, N g{2,... ,p}, (3.43) 
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then G'p+i(/xo, no) 7^ as long as /io ^ {Eq, ■ ■ ■ , E2p+i}- This fact immediately 
follows from (2.16) since Fp(//o.rio) = but -R2p+2(Mo) 7^ by hypothesis. In 
particular, (xj^ino), . . . ,fij^,{na) all meet on the same sheet since 

A>("o) = (Mo,Gp+i(/io,no)), r = l,...,iV (3.44) 

and hence no special divisor can arise in this manner. It remains to study the case 
where two or more projections collide at a branch point, say at {Emo^^) for some 
no e Z. In this case one concludes Fp{z,nQ) = 0[{z — Emo)^) and 

Gp+i{Emo,no) = (3.45) 

using again (2.16) and Fp{Emo,no) = R2p+2iEmo) = 0. Since Gp+i(-,no) is a 
polynomial (of degree p + 1), (3.45) implies Gp+i{z,no) = 0{{z — Emg))- 

Thus, using (2.16) once more, one obtains the contradiction, 

O{{z-Emof) = R2p+2{z) (3.46) 



{z-E^„)i Yl {E„,-Em)+O{z-Emo)]- 



m— 1 
m^mo 

Consequently, at most one fij{n) can hit a branch point at a time and again no 
special divisor arises. Finally, by (3.33), fij{n) never reaches the points Poo+- Hence 
if some fij (n) tend to infinity, they all necessarily converge to Poo_ • Again no special 
divisor can arise in this manner. 

The proof for !?/>(„) is completely analogous (replacing Fp by i?p+i and noticing 
that by (3.33), (f> has no zeros near Poo±)^ thereby completing the proof. CH 

RerriEirk 3.7. For simplicity we assumed a{n) ^ 0, /3(n) 7^ 0, n e Z, in Lemma 
3.6. However, the asymptotic behavior in (3.38) (resp., (3.40)) shows that no special 
divisors can be created at infinity (resp., zero) and hence the results of Lemma 3.6 
extend by continuity to the situation considered in Remark 3.5. In particular, it 
extends to the case where (3{no) = for some no € Z. The case a(no) = for some 
TiQ G Z is more involved and requires more and more case distinctions as is clear 
from Remark 3.5, but the pattern persists. 

Next we turn to the representation of (j), ^, a, and P in terms of the Riemann 
theta function associated with K,p. We freely use the notation established in Ap- 
pendix A, assuming JCp to be nonsingular as in (3.2)-(3.4). To avoid the trivial case 
p = (considered separately in Example 3.13), we assume p e N for the remainder 
of this argument. 

We choose a fixed base point Qo G }Cp\{Poo+,Poo-,Po,+ ,Po,-}, in fact, we will 

choose a branch point for convenience, Qo S BiJCp). Moreover we denote by 

a normal differential of the third kind (cf. (A. 11), (A. 12)) with simple poles at Pi 

and P2 with residues 1 and — 1 , respectively. Explicitly, one computes for _ 

and <jJ^p^ _ the following expressions 

^Po,-,Poo± =^^|^^T^n(^-A±,,)rf^, Po,- = (0,2/0,-) = (0,-5p+i), 

(3.47) 
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where {X±j}j=i,...,p are uniquely determined by the normaHzation 



/ 



'^S-.Pci = 0' i = l'---'P- (3-48) 



The expUcit formula (3.47) then imphes the following asymptotic expansions (using 
the local coordinate C = near Po,± and = 1/z near Poo±)) 

rP 



^ ^pi,p^_ ^Zo {in(c)} +'^o°'^(^o,-,Poo_) +0(C) as P ^ Po,±, (3.49) 
J^^'pI-,p^_ ^=o{-ln(C)}+'^o°°*^^'o,-,Poo_) + 0(C) asP^P,o±, (3.50) 
J^^pi,Poo^ ^=0 {ln(C)} +^o^iPo,-,PooJ + 0(C) as P ^ Po,±, (3.51) 

j^^'pL.P^, ^Zo V^^} +'-o°''(^o.-,PooJ +0(C) as P- Poo^. (3.52) 

Here Qo S S(/C,i) is a fixed base point and we agree to choose the same path of 
integration from Qo to P in all Abelian integrals in this section. 

Lemma 3.8. With oj^" {Pq-, Poo ±) and oJq'' {Pq-, P^c^), a, a' € {+,—}, defined 
as in (3.49)-(3.52) one has 

exp [u)q'~{Po-,Poo±) -t^r"^(-Po,-,-Poo±) 

-a;^-(Po,_,PooJ+a;S'^(Po,-,PooJ] =1- (3.53) 

Proof. Pick Qi,± = {zi,±yi) G /C„\{Poo±} in a neighborhood of Poo± and Q2,± = 
(z2,±y2) G l^n\{PQ.±} in a neighborhood of Po,±. Without loss of generality we 
may assume that Poo+ and Po.+ lie on the same sheet. Then by (3.47), 



^P0,-,Poo_ ~ / ^Po,-,Poc_ ~ / '^Po,-,Poo_ + / '^P0,-,Po< 



Qo •'Qo •'Qo •'Qo 

/■Q2,+ J rQi,+ ij2- 

= ^ = in{z2) - ln(2i) + 2Trik, (3.54) 

^Qo ^ .^Qo ^ 

for some k G Z. On the other hand, by (3.49)-(3.52) one obtains 

/ 

Jc 



'^Po.-.P.o^ " / '^Po.-.Poo. - / '^Po.-.Poo. + , 

Qo -^Qo JQo JQo 

= ln(z2) +ln(lM) +^°'"(^o,-,PooJ -a;o°°+(Po,-,Poo_) -a;o°°-(Po,-,Poo_) 
+ a;S'+(Po,-, Pcx,_ ) + 0{Z2) + 0{\/zi), (3.55) 

and hence the part of (3.53) concerning oj^p^ _ follows. The corresponding 

(3) 

result for ojp^ _ p^^ is proved analogously. □ 
In the following it is convenient to use the abbreviation 

z{P, Q) - Sq„ - (P) + (Dq), P e /Cp, g = {Qi , . . . , Q^} e Sym^ /C^. 

(3.56) 
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Theorem 3.9. Suppose that a, (3 C C satisfy (3.1) and the pth SB system, (2.28). 
In addition, assum,e (3.2)-(3.4) and let P e /Cp \ {f(x)+, -Poo_ , -Po,+, -Po,-} and 
{n,no) € Z^. Then for each n &Z, 'Df^^n) md 2?i/(n) oi'^e nonspecial. Moreover,^ 

MP, n, no) = C{n, no)^^^^p-^ exp [{n - no) a;i,j,_,p„ J , (3.58) 
V'2(-P, n, no) = C{n)C{n, uq) 



where 

C(n) = (-l)"-"''exp [(n-no)(a;S'-(Po,-,P<x>_) -a;o°°+(Po,-,Poo_))] 

Q;(no) " 0(z{Po.+,ii(na))) 

^(z(Poo+,A(no))) ,,,,, 
C(n, no) = exp [ - (n - no)a;o +(Po,-, Poo+)J ^(^(p^ "A(n))) " ^ ^ ^ 

T/ie ^6e/ map linearizes the auxiliary divisors in the sense that 

= ^QoC^Kno)) + Ap„__iPoo+)in - no), (3.62) 

^Qo(%«)) =^iQo(%no)) +^Po,-(^oo+)('^-'^o)- (3.63) 

Fi'lially, a, (3 are of the form 

a{n) = a(no) (-!)"-"« exp [ - (n - no)(c^o°'"(^o,-, Poo_) - a;o°°+(Po,_, PooJ)] 
^(^(Po,+,^(no)))0(z(Po,+,A(n))) 



e(^(Po,+, A(no)))^(^(Po,+, H"))) ' 
/?(n) = /3(no)(-l)"-"'' exp [(n - no)(a;S'-(Po,-, PooJ - a;o°°+(^'o,-,Poo_))] 
^ g(^(Poo+ , Mno)))g(^(Poc^ , u{n))) 

e{z{Po.^,P{m)))6{z{Po.^,fi{n))y 

a{n)p{n) = exp [lu^^ (Po,-, PooJ - co°'+{Po,-,PooJ] 
^ 0{z{Po,+,Kn)mz{P^^,v{n))) 
0{z{Po.,+,iL{n)))e{z{P^^,ii{n))y 



(3.64) 



(3.65) 



(3.66) 



Proof. While equation (3.62) is clear from (3.36), equation (3.63) follows by com- 
bining (3.15) and (3.37). By Lemma 3.6, "D^ and "D^ are nonspecial. By (3.15) and 

Theorem A. 3, (j){P, n) exp ( — ^^p} _,p^ ) must be of the type 



^(P,n)exp ( - / .(^)_,,^ ] = Cin) ^^ (3.67) 



^To avoid multi-valued expressions in formulas such as (3.57)-(3.59), etc., we always agree to 
choose the same path of integration connecting Qo and P. 
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for some constant C{n). A comparison of (3.67) and the asymptotic relations (3.33) 
and (3.34) then yields the following expressions for a and /3: 

= -Ce"o°°-(^o,-.P^_) ^fe(^oo_,^)) (3.68) 

%(-P<x>_,M)) 

Similarly one obtains 

^+ ^+.c.r+fPn.-,P„. )^(^(^oo+,^+)) 



Here we used 

aQo(^A+) =^Qo(^a) +^Po,-(^oo+), aQo(^^+) =^Qo(%) +4Po,-(^ooJ, 

(3.70) 

(3.56), and relations of the type 

z{Poo+,[i^) = z{Poo.,K)=ziPo,-,ti) = ziPo,+,E^), (3.71) 

z{Poo+,E^) = z{Po,-,0), 0(Po,+ ,A+) =^(^'<x,_,A)- (3.72) 

Thus, one concludes 

C(n + l) = -exp[a;°'-(Po,-,fcx>_)-w^+(Po,-,i'oo_)]C(n), neZ (3.73) 

and 

C7(n + l) = -exp[a;^-(Po,-,Poo_)-wS'^(fo,-,Pcx=_)]C(n), n e Z, (3.74) 

which is consistent with (3.53). The first-order difference equation (3.73) then 
implies 

C{n) = (-l)(-»o) exp [(n - no){to°'- {Po,-,Poo. ) - co^-" (^o,-, Poo_ ))] CM, 

n, no G Z. (3.75) 

Thus one infers (3.64) and (3.65). Moreover, (3.75) and taking n = no in the first 
line in (3.68) yield (3.60). Dividing the first line in (3.69) by the first line in (3.68) 
then proves (3.66). 

By (3.36) and Theorem A. 3, ipi{P, n, no) must be of the type (3.58). A compar- 
ison of (3.18), (3.33), and (3.58) as P ^ Poo+ (C = ^/z) then yields 

Vi(P,n,no) r°-"(l + 0(C)) (3.76) 

and 

, . ^, , ^(l(Poo+,A(n))) 

ipi{P,n,no) = C{n,no)—7^ — v^- 

' C^o ' a-Poo+,M^^o))) 



0{ziP^+ 


,A+)) 


0{z{Po,- 


,0) 


OiziPo,- 


,A)) 


0{z{Po,- 




0{z{Po,- 
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X exp [{n - noH°°+(Po,-,Poo+)]C"''""(l + 0(C)) (3.77) 
proving (3.61). Equation (3.59) is clear from (3.21), (3.57), and (3.58). □ 

Remark 3.10. (z) By (3.62), (3.63), the arguments of all theta functions in (3.57)- 
(3.59) (3.61), and (3.64)-(3.66) arc linear with respect to n. 
{ii) Using relations of the type (3.71), (3.72) and 

«Qo (^^) = (^Qo (^a) + Ap„,_ (Poo^ ), (3.78) 

one can rewrite formulas (3.57)-(3.66) in terms of p, (or only. 
(Hi) For simplicity we assumed a{n) ^ 0, /3(n) 7^ 0, n e Z, in Theorem 3.9. Since 
by (3.33) and (3.34) no jlj and no i)£ hits Po,+ or P00+, the expressions (3.64) and 
(3.65) for a and f3 are consistent with this assumption. 

(iv) Generally, a and (3 will not be quasi-periodic with respect to n e Z. Only under 
certain restrictions on the distribution of {£'m}m=o 2p+i! such as the (de)focusing 
cases discussed in Corollary 3.11 next, one can expect to uniformly bound the 
exponential terms in (3.64) and (3.65) and prove quasi-periodicity of a and (3. 

The special defocusing and focusing cases are briefly considered next. 

Corollary 3.11. Suppose that a, /3 C C satisfy (3.1) and the pth SB system (2.28) 
and assume (3.2) -(3.4). Moreover, assume either the defocusing case, where /3(n) = 
a{n), or the focusing case, where (3{n) = —a{n), n £Z. In either case, a is quasi- 
periodic with respect to n gZ. 

Proof. We start by noting that the ratio of theta functions in (3.64) and (3.65) 
is bounded as n varies in Z since by (3.15) (see also (3.33) and (3.34)) Po,+ is 
never hit by any i>e{n) and P00+ is never hit by any fij{n). Thus, a (and of course 
/?) is quasi-periodic if and only if the exponential term in (3.64) is bounded (i.e., 
unimodular). Assume the defocusing case (3 = a. Then, writing 

a{n) = 6(n)e"^ /3{n) = 6(^^)e-"^ n e Z, b,b € £°°{Z) (3.79) 
(cf. (3.64), (3.65)), 13 = a implies 

(3{n) = ^(rj)e-"i^«(^)-'"i'"('^) = a{n) = ^e"^'^('=)-'"^"'(^) (3.80) 

and hence Re(c) = 0. The analogous argument applies in the focusing case. □ 

Remark 3.12. (i) The additional (de)focusing assumption f3 = ±75 in Corol- 
lary 3.11, implies strong restrictions on the possible location of the branch points 
{Em, 0), m = 0, . . . , 2p + 1. In particular, in analogy to the Ablowitz-Ladik model 
discussed in [34], one expects all {Erm 0) to occur in pairs which arc reflection sym- 
metric with respect to the unit circle T in C. In the defocusing case, (3 = a with 
|a(n)| < 1, n e Z, all branch points are seen to lie on T as discussed in [22] and 
[38]. For |q;| > 1 one expects them to bifurcate off the unit circle T. 
{ii) In analogy to the defocusing case of the nonlinear Schrodinger equation (cf. 
[27, Ch. 3]), the isospectral manifold of algebro-geometric solutions of (1.3) can 
be identified with a. (p + l)-dimensional real torus T^'"'""'^ as discussed in detail in 
[38, Ch. 11]. This isospectral torus is of dimension p+1 (rather than p, given the 
p divisors fij{no), j = 1,. . . ,p) due to the additional scaling invariance discussed 
in (2.30), (2.31) involving an arbitary constant multiple of absolute value equal to 
one. 

{Hi) By Remark 3.7, no special divisors arise if /3(ni) = ±a{ni) = for some ni e Z 
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and hence Corollary 3.11 extends to this case as long as (3{no) = ±Q;(no) 7^ in 
(3.65). 

{iv) In the special defocusing case P = a, with \a{n)\ < 1, n G Z, Corollary 3.11 
recovers the original result of Geronimo and Johnson [22] that a is quasi-periodic 
without the use of Fay's generalized Jacobi variety, double covers, etc. 

Finally, we briefly consider the case p = excluded in Theorem 3.9. 

Example 3.13. Letp = 0, P = {z,y) e /Co\{Po,+ , -Po,-, -P00+, -Poo_}, and (n,no) e 
Z2. Then, 

/Co : M^, y) = y^- R2{z) = y^-{z- Eo){z - Ei) = 0, 

Eo,E^ e C\{0}, So ^ El, gl = EqEi, gi = y(Po,+), ci = -{Eo + Ei)/2, 
a{n) = a(no)(-gi)"-"», P{n) = /?(no)(-5i)"°-", 

s-SBo(a, /3) = f L-V/m^) = 0' = [1 - (ciM)]/2, 



2(/3-+ffi/3) 
y + z- 2a+p + ci -2pz 



-2a+ y - ^ + 2a+/? - ci ' 

One verifies that Eq ^ Ei is equivalent to a/3 € C\{0, 1}. For a Borg-type 
theorem related to this example in the special defocusing case (5 = a with \a{n)\ < 
1, n e Z, we refer to [31]. 

Appendix A. Hyperelliptic Curves and Their Theta Functions 

We give a brief summary of some of the fundamental properties and notations 
needed from the theory of hyperelliptic curves. More details can be found in some of 
the standard textbooks [19] and [35] , as well as monographs dedicated to integrable 
systems such as [14], Ch. 2, [27], App. A, B. 

Fix £( € N. The hyperelliptic curve fCg of genus g used in Section 3 is defined by 

2g+l 

Kg : Tg{z, y) = y^- R2g+2{z) = 0, R2g+2{z) =Jl{z- E^), (A.l) 

{-E„i}m=o,...,2g+i C C, S„ 7^ for m 7^ to', m,m' = 0, ... ,25 + 1. (A. 2) 

The curve (A. 2) is compactified by adding the points P00+ and Poo- , P00+ Poo- > 
at infinity. One then introduces an appropriate set oi g + 1 nonintersecting cuts Cj 
joining -5^0) and Ej^'Q)- We denote 

C= U Cj, CjnCk = 9, ji^k. (A.3) 

iG{l,...,s+l} 

Define the cut plane 11 = C\C, and introduce the holomorphic function 

/2g+l \ 1/2 

i?2s+2( • Y'^ : n ^ C, z^{^{z-Em)\ (A.4) 

\m=0 / 

on n with an appropriate choice of the square root branch in (A.4). Define 

Mg = {{z,aR2g+2{zf'^)\z&C, ere {±l}}U{Poo+,Poo_} (A.5) 
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by extending i?2g+2( • )^/^ to C. The hyperelliptic curve K,g is then the set Mg with 
its natural complex structure obtained upon gluing the two sheets of M.g crosswise 
along the cuts. The set of branch points B{JCg) of Kg is given by 

BiKg) = 0)}„=o,...,2s+i (A.6) 

and finite points P on ICg are denoted by P = {z,y), where y{P) denotes the 
meromorphic function on ICg satisfying J^g{z,y) = — R2g+2{z) = 0. Local coor- 
dinates near Pq = {zq, yo) e Kg\{B{Kg) U {Pcx)+ , P<x>- }) are given by Cpq = z - zq, 
near Poo-t by Qp^^ =\lz, and near branch points (i?m„,0) G B{JCg) by C{E^g,o) = 
{z — Emoy^"^- The Riemann surface ICg defined in this manner has topological genus 
9- 

One verifies that dz/y is a holomorphic differential on ICg with zeros of order 
5 — 1 at Pooj. and that 

z^~^dz , , , . „, 

nj = ——, j = l,...,g (A.7) 

form a basis for the space of holomorphic differentials on ICg. Introducing the 
invertible matrix C in C^, 



C — {Cj,k)j.k=l,...,g, Cj^k — I Vj 



(A.8) 



c(fc) = (ci(fc),...,c„(A:)), Cj{k) = C^l, j,k = l,...,g, 
the corresponding basis of normalized holomorphic differentials Uj, j = 1, . . . ,g on 



ICg is given by 

g „ 
Wj ='^Cj{e)rie, / ojj=Sj,k, j,k = !,...,( 



(A.9) 



Here {aj, &:/"}j=i,...,g is a homology basis for ICg with intersection matrix of the cycles 
satisfying 

cij obk ^ Sj^k, aj o afe = 0, bj obk = 0, j,k = 1, . . . ,g. (A.IO) 

Associated with the homology basis {aj,bj}j=i^...^g we also recall the canonical 
dissection of ICg along its cycles yielding the simply connected interior ICg of the 
fundamental polygon dICg given by dICg = aibia^^bi^a2b2a2^b2^ ■ ■ ■ aj^b~^. Let 
M{ICg) and M^{ICg) denote the set of meromorphic functions (0- forms) and mero- 
morphic differentials (1-forms) on ICg. The residue of a meromorphic differential 
ly G Ai^ilCg) at a point Q E ICg is defined by resQ(i^) = J^^ i', where -jq is a 
counterclockwise oriented, smooth, simple, closed contour encircling Q but no other 
pole of v. Holomorphic differentials are also called Abelian differentials of the first 
kind. Abelian difi^erentials of the second kind, w^^^ € M.^{ICg), are characterized by 
the property that all their residues vanish. Any meromorphic differential lo^^^ on 
ICg not of the first or second kind is said to be of the third kind. A differential of the 
third kind w^^^ G M^{ICg) is usually normalized by the vanishing of its a-periods, 
that is, 



/ a;(3)=0, j = l, 



(A.11) 
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A normal differential of the third kind associated with two points Pi, P2 € 

^5, -Pi 7^ P2 by definition has simple poles at Pj with residues (—1)-'+^, j = 1,2 
and vanishing a-periods. If Wp_Q is a normal differential of the third kind associated 
with P, Q G )Cg, holomorphic on )Cg\{P, Q}, then 

where the path from Q to P lies in ICg (i.e., does not touch any of the cycles aj, 

We shall always assume (without loss of generality) that all poles of differentials 
of the second and third kind on ICg lie on ICg (i.e., not on dlCn)- 
Define the matrix r = (rj,^)j,£=i,...,g by 

Tj,i= <^j' J,^=l,---,5- (A.13) 

J be 

Then Im(T) > and Tj^e = t^j, j,i = 1, ... ,5. Associated with r one introduces 
the period lattice 

Lg = {zeC^ \z = m + nT, m,n&'^^} (A.14) 

and the Riemann theta function associated with /Cg and the given homology basis 
{aj, 6j}j=i,...,g, 

e{z) = ^ exp {2Tri{n,z) + m{n,nT)), z e C, (A.15) 

where {u,v) = ^vJ = Sf=iW^j denotes the scalar product in C^. It has the 
fundamental properties 

e{zi, . . . , Zj-i,-Zj, Zj+i, . . . , Zg) =6{z), (A.16) 

^(^ + m + nT) = exp ( - 27ri(n,^) - 7ri(n,nT))^(^), m,nGlI>. (A.17) 

Next, fix a base point Qo G ^ff\{-fo,±, -Poo-t}, denote by J{lCg) = C^/Lg the 
Jacobi variety of ICg, and define the Abel map Aq^ by 

Aq, : tCg ^ JitCg), Aq^{P)=( f . . . , / (mod Lg), P G )Cg. 

\JQo JQo / 

(A.18) 

Similarly, we introduce 

Oq^ : Div(/Cg) ^ J{lCg), V ^ {V) = ^ V{P)Aq^ (P), (A.19) 

where Div(/Cg) denotes the set of divisors on Kg. Here I>: /Cg ^ Z is called a 
divisor on /Cg if T>{P) 7^ for only finitely many P € ICg. (In the main body of this 
paper we will choose Qo to be one of the branch points, i.e., Qo S S(/Cg), and for 
simplicity we will always choose the same path of integration from Qo to P in all 
Abelian integrals.) 

In connection with divisors on /Cg we shall employ the following (additive) no- 
tation, 

T^QoQ = T^Qo + Vq^Vq,+---+ Vq^ , (A.20) 

Q = e Sym™/Cg, Qo e /Cg, m e N, 
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where for any Q e /Cg, 

and Sym" /Cg denotes the nth symmetric product of ICg. In particular, Sym"'/Cg 
can be identified with the set of nonnegative divisors < V G Div(/Cg) of degree 
m. 

For / G M{ICg)\{0}, Lo e Al^(/Cg)\{0} the divisors of / and uj are denoted 
by (/) and (w), respectively. Two divisors V, £ Div(/Cg) are called equivalent, 
denoted hyV^S, if and only if V-S = (/) for some / e M{JCg)\{0}. The divisor 
class [V] of V is then given by [V] ^ {£ e Div(/Cg) \ £ ^V}. We recall that 

deg((/)) = 0, deg((c.)) = 2(5 - 1), / e M(/Cg)\{0}, uj e Mi(/Cg)\{0}, (A.22) 

where the degree deg(I?) of T> is given by deg(I?) = J2peK '^i-^)- if) called a 
principal divisor. 

Introducing the complex linear spaces 

>C(D) = {/eM(/Cg) |/ = 0or (/)>P}, r{V)=dimC{V), (A.23) 

C^{V) = {ojGM^{ICg)\(j = Oor {oj)>V}, = dimC^{'D) (A.24) 

with the index of speciality of V, one infers that dcg(X>), r{'D), and i{V) only 
depend on the divisor class [D] of T>. Moreover, we recall the following fundamental 
facts. 

Theorem A.l. Let V e Div(/Cg), w e M^iKg) \ {0}. Then 

liV) = rip - (w)), g e No. (A.25) 
The Riemann-Roch theorem reads 

r{-V) =degiV)+i{V)- g+1, .9 G Nq. (A.26) 
By Abel's theorem, V e Div(/Cg), 5 G N, is principal if and only if 

deg(I?) = and oq^ {V) = 0. (A.27) 

Finally, assume g G N. Then Oq^ : Div(/Cg) ^ J{^g) is surjective [Jacobi's 
inversion theorem). 

Theorem A.2. Let Vq e Sym^/Cg, Q = {Qi, . . .,Qg}. Then 1 < ^(Dq) = s < 
(7/2 i/ anrf only if there are s pairs of the type {P, P*) € {Qi, ■ ■ ■ , Qg} {this includes, 
of course, branch points for which P ~ P*). 

Denote by Sg^ = (S 

Qo,n ■ ■ • I '^Qo.g) vector of Riemann constants, 
^Qoj = ^(1 + Tjj) - / MP) / a;,-, j = l,...,g. (A.28) 

Theorem A. 3. Let Q = {Qi, . . . , Qg} G Sym^ /Cg and assume Vq to be nonspecial, 
that is, HVq) = 0. Then 

0{^Q, - AQoiP) + aQoil^o)) = Oif and only if P G {Qi, . . . , Qg). (A.29) 
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